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SUMMARY 


A theory is developed for aircraft control obtained from the propeller 
forces and moments generated by blade angle variation during a blade revolution. 
The propeller blade is pitched harmonically one cycle per propeller revolution. 

As a result, the propeller off-set thrust and power produce forces and moments 
on the axis which provide control for the vehicle. This is termed cyclic-control. 
Using quasi-steady principles, and a power series representation of an arbitrary 
function for cyclic-blade angle, theory is developed which leads to exact solu- 
tions in terms of derivatives of steady-state thrust and power with respect to 
blade angle. An alternative solution, when the cyclic-blade angle is sinusoidal, 
is in terms of Bessel functions. This solution is mathematically valid for 
cyclic-blade angles to 90 degrees. Steady-state propeller theories and experi- 
ments are well documented and available in the literature, and in addition, an 
approximate method is given in the present report. Needed are the steady-state 
first and higher order derivatives of thrust and power with respect to steady- 
state blade angle. An estimate of non-steady azimuth angle change or lag is 
presented. Correlations of cyclic-control theory predictions are made with 
experiment. 

For vehicle control it is desirable to isolate forces and moments. In 
this aspect it is shown how the counterrotating propeller is unique. By 
applying a cyclic-blade angle at the proper azimuth angle on each rotor, a 
force can be generated in any direction without an accompanying moment on 
the counterrotating propeller axis, or alternatively, a moment without force. 

In the hovering mode, cyclic control provides propulsion in any direction to 
the dual rotor without rotor tilt or moments. Also, when in the propeller mode 
at the tail of an airship or submarine, this provides isolated or combined con- 
trol forces and moments with no airship movement. It is shown how three or 
more single rotation propellers, in line as on a wing, can also isolate control 
forces and moments by cyclic-control. 


INTRODUCTION 

In the technical evolution of broadening the high-low speed range capa- 
bility of aircraft, the propulsive unit, as speed decreases, must increasingly 
augment both lift and control. Ideally, both would be contained in a single 
and optimized unit. By making the blades cyclic controllable, the propeller 
or rotor becomes a propulsive unit in which propulsion, lift, and control are 
combined, as in some helicopter designs. The term cyclic control refers to 
harmonically changing the blade angle from Bo + y to Bo - Y to e 0 + Y degrees 
during each propeller revolution, where y and the azimuth angle, v 0 , at which 
y is added, are both controllable. A list of the types of aircraft or the 
types of operations for cyclic control application would include tilt- 
wing, tilt-rotor, vectored slipstream by normal or side forces, short take-off 
and landing, loiter, lifting-crane, airship or submarine, control of slipstream 
deflection for best cruise, and air-cushion. In some of these applications 
only a limited number of degrees of freedom is of concern and the term mono- 
cyclic propeller has been used. Cyclic-control theory does not become more 
complicated for more degrees of freedom and can be derived for generalized 
cyclic control valid for six degrees of freedom. 



Cyclic-control experimental data is presented in reference 1 on work 
done by de Decker. This experimental data is with zero free stream velocity 
and includes thrust and power, without and with cyclic blade-angle, and side 
force and pitching moment of the propeller with cyclic blade angle for various 
blade angle settings. Basic cyclic-control analysis is developed in reference 
2. Recent counterrotating- rotor cyclic-control analysis for the helicopter, 
based on helicopter theory assumptions of low disc loading and small blade and 
inflow angles, was developed by Drs. Wadi a and Fairchild and presented (ref. 3) 
in a paper. The analysis of reference 2 is partly based on an analogy of cyclic- 
control, with inclined propeller methods. The theory of a propeller with cyclic 
change of blade angle is analogous to the theory of the inclined propeller or 
the propeller in yaw developed in references 4 and 5. The primary difference 
is that the cyclic-control blade does not sense a dynamic pressure change, 
whereas, for the inclined propeller, changes in dynamic pressure occur during 
a cycle. Helicopter theory for cyclic control is simplified to rotors of small 
solidity and small blade and inflow angles, and is generally inadequate for 
highly loaded propellers of large solidity. The theory of reference 2 contains 
the concept that cyclic control on a counterrotating propeller can isolate all 
control forces and moments. Details of this concept were presented by DeYoung 
for VLM Corporation to an Army evaluation team USAAMRDL (Ames Research Center), 
in an unpublished paper, "The Counterrotating Rotor with Cyclic-Control Permits 
Isolation of Forces and Moments," February 17, 1972, Fort Worth, Texas. 

The principal objective of the present study is to present a generalized 
theory relating cyclic-control forces, moments, power and thrust changes to the 
steady-state characteristics of the same propeller for arbitrary propeller 
solidity or disc loading; to derive counterrotating propeller cyclic-control 
analysis in terms of that of the single propeller solution; and to investigate 
the effectiveness of propeller or rotor configurations or combinations for 
isolating cyclic-control forces and moments. The present study starts with a 
refinement, clarification, and extension of the basic theory developed in 
reference 2. 


SYMBOLS* 


B 

b 


number of blades 


propeller blade chord (b gv = 
pitching moment coefficient 



b dp) 


M 

p n^ D^ 

a 


* All forces and moments are acting at propeller hub. See figure 1 for 
positive directions. 


2 


normal -force coefficient 
yawing moment coefficient 


power coefficient (2ttCq), 
torque coefficient (Cp/2iT 
thrust coefficient 


side-force coefficient 


N 

Nyaw 
Pa " 2 ° 5 




), 




„2 n* 

p a " D 



D 


4 


propeller diameter 

figure of merit (static conditions. 



advance ratio (V/nD) 

advance ratio at which thrust is zero 

advance ratio at which power is zero 

Bessel function of the first kind (eq. (44)) 

fractional change of power radial position with blade angle 
(eq. 57)) 

fractional change of thrust radial position with blade angl 
(eq. (53)) 

propeller rotational speed, revolutions per second 


power input to propeller 


2 

freestream dynamic pressure (p„V /2) 

a 

propeller radius 

radius distance from propeller axis 
propeller thrust 
freestream velocity 
blade angle at p = 3/4, (g - e) 

blade angle at p = 3/4, measured, from zero lift (e + e ) 

blade angle at p = 3/4 for which p tang(p) is most constant along 
radius, measured from zero lift {eq. (54)) 

maximum cyclic angle, positive at v = it to give positive pitching 
moment 

blade angle between zero lift chord and geometric chord 
section Co /2tt 

dimensionless radial distance (r/R) 
density of air 

thrust off-set for pitching (2 C/C T ) 

m i av 

radius center-of-power position (r p /R) 

Pp at blade angle g Q = e-j (eq. (57)) 

radius center-of-thrust position (r-p/R) 

p t at blade angle 6 0 = e-j (eq. (53)) 

48 b av 

effective solidity ( 3 ^ D ) 


T 


blade tip correction factor (eq. (66)) 


<j> propeller inflow angle 

y azimuth angle, positive clockwise 

azimuth angle for maximum cyclic-blade angle (eq. (1)) 

Subscripts 
av average 

o indicates the steady- state value 

Superscripts 

first, second, etc., derivatives with respect to or fj 
* denotes linear with advance ratio from Jqj 



SOLUTION FOR CYCLIC-CONTROL PROPELLER CHARACTERISTICS 


Force and Moment Coefficients in Terms 
of Azimuth Angle Integrals 

For cyclic-control of propeller forces and moments a cyclic variation of 
blade angle is superimposed on the steady-state blade angle during one 
revolution of the propeller. When this variation is continually repeated, the 
thrust and torque will also vary during the cycle and lead to thrust off-set, 
side forces, and resulting moment changes. The azimuth angle, ¥, is defined 
as shown in figure 1, also the positive directions of force and moment 
coefficients are shown. Let 

for clockwise rotation; h , q = y ^ + 4» ax + y con 

for counterclockwise rotation; v = -y . + v + f „ 

o ph ax con 

where is the azimuth angle at which the blade thrust and torque are minimum, 

fph is 0 the phase angle that the thrust and torque lag the blade angle, Tax N 
is the effective phase angle change due to blade pitching velocity (appendix A), 
and T CO n is the azimuth angle that can be chosen to govern or control the 
position of minimum blade thrust. The average thrust and average power 
coefficients are 

2tt 

c T (T“T 0 )d’if (2) 




C Pav 



2t r 

c p(T-T 0 )d¥ 


(3) 


where Cj(v-v 0 ) and CplT-^) represent functional relations of variations with 
(f-^o), and are abritrary functions. 

The hub pitching moment is the product of the thrust and the z component 
of radial distance of center of thrust, then 


^m " 4 tt 


p t(t-t o ) c t(t-4'o) co s v dy 
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The hub yawing moment is the product of the thrust and the y component of 
radial distance of center of thrust, then 


f 2ir 

P T (M 0 )C T (M 0 )sin * dr 


= C tan ¥_ 
m o 


(5) 


where the second equality is obtained by noting that the functional variation 
of pj and Cj is antisymmetric with (¥-¥ 0 ) and that 

sin ¥ = cos ¥ q sin(¥-¥ Q ) + sin ¥ Q cos(¥-¥ Q ) 


then the first term in this sin ¥ expansion will be zero in the integration 
from ¥ = 0 to 2-rr. Similarly, the expansion of cos ¥ in equation (4) leads to 
cos ¥ 0 cos(¥-¥ 0 ) and this multiplied by tan ¥ 0 results in C n = C m tan ¥ 0 . 

The side force is the integration of the change of torque normal to the 
xz-plane, divided by the radial center of torque or power, then 


C 


Y 



Pp(V-¥ 0 ) 


COS ¥ d¥ = - 



r 2 'c p (»-* 0 ) 

p p(' F “ ¥ 0 ) 


cos 


¥ d¥ 


( 6 ) 


The normal force is the integration of the change of torque normal to the 
xy-plane, divided by the radial center of power, then 


C 


N 




Cp(t-¥ 0 ) 

Pp(f-Y 0 ) 


sin ¥ d¥ = 


-Cy tan ¥ q 


(7) 


where as in equation (5) the integrand factor of sin ¥ is antisymmetric. 
The thrust off-set for pitching is 


p m = 

m 


M 


2 C 


m 


R T 


av 


"Tav 


(8) 
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The thrust off-set for yawing is 


p n = 


f T , 

i yaw 


2 C 


n _ tan 
— = m o 


"Tav 


The figure of merit (static conditions) is 

r 3/2 
l T 



(9) 


00 ) 


The signs of equations (4) through (7) are for clockwise propeller 
rotations, however, the final signs and magnitudes of these coefficients will 
depend on the azimuth angle, Yq. For counterclockwise rotation, equations (4), 
(5), and (7) remain the same, only the sign of equation (6) for Cy changes, but 
the final sign and magnitude will depend on f 0 . With the azimuth control 
angle, y con , equal to zero, and T ax small, then the counterclockwise rotation 
has C m and C^j with the same sign, and Cy and C n with opposite sign to that at 
clockwise rotation, with magnitude the same. 


Cyclic-Control Solutions by Power Series 

All the functional factors in the integrands of equations (2) through (7) 
are functions of the blade angle which in steady-state operation would be 
constant with respect to azimuth angle and equations (4) through (7) would be 
zero. In cyclic control operation the blade angle is made to vary during the 
revolution. Let the blade angle be represented by the function 

B(v-h» 0 ) = + yf(f-4' 0 ) (11) 

where y is the maximum value of the cyclic angle, and f(4'-y 0 ) represents the 
variation antisymmetrically of the cyclic angle added to the blade angle during 
a revolution of the blade (e.g. this function is simply -cos(h’-T 0 ) for 
sinusoidal variation of blade angle). Then the thrust coefficient for a given 
blade rotational speed and advance ratio, varies as 


Ct^-Vq) = c t [b 0 + Y f(y-H' 0 )] (12) 

where the right side denotes that Cy is a function of the term inside the 
bracket. Equation (12) can be expanded into a power series in terms of y, thus 


8 


C T = C T 


♦ ^i) y + — ?!fi) y 2 + _L !!!ib + 

Y=0 9y / y=0 2! 8Y ' y=0 31 9Y l Y =0 


The coefficients of the y series are obtained by taking partial derivatives 
of equation (12), and letting y be zero. Then 


C T = c t C 3 0 + = C T 


9C t 9C t Y=0 9Ct 

_ L = ! f(y-v) = 5. f (v - Y ) 

3y 3[6 0 + 9B 0 


9 C T 9 C T ~ Y=0 9 Ct 0 

-p l = 1 7 [f(*-0] = -4^ [f(*-0] 2 

3Y 2 3 [B q + Yf('i'-T 0 )] 2 ° 3B 0 2 ° 


Then equation (12) can be written as 


3Cy Y 2 3 2 Ct 0 

r f(*-f 0 ) + y [f(T-T n )] Z 

ae„ 0 2! nj 0 


y 

— f (V 

9f3 o 


3 

3 ~ 

Y 

9 C T 0 

3! 

«o 3 

Y n 

3 " c t„ 


n! 36. 


r + 


when Cy Q is the steady-state thrust coefficient (y = 0) which is a function 
of 3 0 . Thus in equation (14) the cyclic variation of thrust is expressed 
simply in terms of steady-state thrust derivatives with respect to blade 
angle. 
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A similar expansion gives the power coefficient as 


c p(f-f 0 ) 


aCp y 2 a 2 Cp 2 

Cp + Y f (¥-¥_) + — 2^ C f ( 4 '-' 1 'o)3 + 

P o 96 0 0 2! 98o 2 


3 ,3- 


Y <3 J 

-2 [f(T-Tj] + 

3! 3B 0 3 


Y n 3 n Cp n 

+ J^[f(*-0] + 

n! 8S 0 


An expansion for the product of radius center-of-thrust and thrust, 
results in 


3 ( p To C To) 

Pt^-O C t (¥-¥ ) = P T C T + Y — H'V-'Vq) + 

101 0 0 3B„ 


Y 2 ^(pTqCtJ 2 

+ 

2! ae 0 2 


Y n 9n ( p T 0 C T 0 ) r , " 

+ 1 — [f(f-%)] + . . . 

n! 3* 0 


An expansion of the ratio of power to radius center-of-power gives 


C P. 

9 (sr 2 -) 


.2, C Po, 


C p (T-T ) Cp 9( »p 1 Y 2 * ( »p' 2 

— — = ~ + Y - 2 - f(T-T 0 ) + j 2 - [Hl'-'i'J] * 

~ f \u \u \ 'io O \ 'S o 


P< M 0 ) P 0 


2 ! 36 , 


t" 8 " ( %> 

flflHj] + 

n! 36 0 " 
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Inserting equations (14) through (17) into the corresponding integrands 
of equations (2) through (7), and noting that f('P-'F 0 ) is antisymmetric so 
many of the integrals are zero, leads to the solutions 
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I 


where the I- integrals are defined as 


I 


1 r 2lT 

* X 


Pn - I [f<*-*o>3 n d " 
* J o 


(24) 


Mn 


r<- T* 


[ f ( ) ] cos v dY 


(25) 


Nn 


/-2tt 

J r\ 


[f(v-^ 0 )] n sin v dv = I Mn tan v Q 


n=odd 


(26) 


In equations (18) through (23), the single prime indicates the first partial 
derivative with respect to S 0 , the double prime indicates the second partial 
derivative with respect to B 0 , etc.. Equations (18) through (23) are 
mathematically exact series but remain generalized in that an arbitrary cyclic j 
antisymmetric function of y, represented by f(y-'t'o), can be inserted through J 

the I-integrals of equations (24) through (26). The factors of y n in the series ; 
are derivatives with respect to steady-state blade angle of thrust and power i 

coefficients, or of products or ratios of thrust and power with blade radius 
center of thrust or power. 

Cyclic function of cosine to m-power . - The cyclic function giving a i 

sinusoidal variation and distorted sine variations of blade angle are of 
immediate interest. Let j 

f (v-v 0 ) = cos m ('F-^ o ) (27) 


where m denotes odd integer powers or odd integer roots. The integrals of 
equations (24) through (26), with equation (27) are 


1 



IT 



cos a (y-'t' Q ) dy = 


0 , a = odd 

4 r(^^-) , a = even 

(a+1 )i/iT r ^ a+2 ^ 


(28) 
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‘Mn 


, r2, 

7T J 


cos (^-^ 0 )cos^ dm = 


2 r(^) 

cos v, a = odd 

ViT „,a + 3, 0 


r(^) 


( 29 ) 


0 


a = even 


Nn 


■rf 


2tt 


cos (y-^sirw dY 


' 2 r(^) 

— sin <f, a = odd 

vtt r,a+3v 0 

l 2 ; 


a = even 


(30) 


where a=nm, r is the Gamma function, and a = odd, means odd integers or 
odd fractions, 1/3, 1/5, etc.. Using various values of m in equation (27), 
the effect of the function f^-^) can be assessed by evaluating the integral 
factors in equations (24) through (26). With equations (28) and (29), the 
integral factors for the first two terms in the series, are 


m 

cyclic-curve 

description 

cos m (^ 0 ) 

factors for 
Cj or Cp 

eq. (28) 

factors for 

C m or C Y 
eq. (29) 

l P2 

X P4 

T M1 

X M3 

1 

i 

sinusoidal 

1.00 

0.75 

1 .00 cos ¥ 

0.75 cos 

1/3 

bulged-out 

1.43 

1.16 

CT> 

O 

O 

CO 

o 

1 .00 cos 4 ' 0 

1/5 

bulged-out 

1.60 

1.36 

1 .20 cos v 

0 

1 .09 cos 'i* 

0 

step-function 

2.00 

2.00 

1 .27 cos ¥ 

1 .27 cos f 

3 

reflex 

0.63 

0.45 

0.75 cos v 

0.49 cos 




_i 

0 

0 


The integral factors for C n or -C^ determined from equation (30), are the same 
as those for C,^ or Cy» except sin Y 0 replaces cos t 0 , that is, the tan f 0 
relation of equations (21) and (23) holds true. As an example, the integral 
factors inserted into equations (18) through (23), for the sinusoidal curve, 
m = 1 , appear as 


13 




r r _ 1 r 11 2 . .75 r ,v 4 

C Tav ' C To 4 C T Y 48 C T Y 


r r _ 1 r " 2 . .75 r ,v 4 

C Pav " C P Q 4 C P Y 48 C P Y 


C m "4 WV y cos 24 y3 cos 


C tan v n 
m o 


> (31) 


Cy - 2 tt (^p^Pp) Y ^os ~ 1 2 tt (^p^Pp^ Y COS 


”^Y ^ an 


Examination of the Imi factors in the table shows that the step-function 
cycle generates the largest moments and forces, about 27 percent greater than 
that of the sinusoidal cycle. However, as shown by the Ip2 values, the extra 
power required is double that of the sinusoidal cycle. In the step-function 
cycle the maximum blade-pitching acceleration becomes infinite. 


Cyclic function for most moment with least extra power, also, with minimum 
blade-pitching acceleration . - For cyclic control an objective is to get the 
most control of forces and moments with the least extra power requirement. A 
relationship between moment or force and extra power can readily be obtained 
from equations (19), (20), and (22). Up to moderate values of the y-angle, the 
second integral terms in equations (18) through (23) are negligible, then 
equation (19) for extra power can be written as 


Solving for y 


Y= 2 


C P 

0 

C P 

4 

Y 


( Cp av 

1 

O 

"O 

o 

W 

! P2 / 


Id9 
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then with equations (20) and (22) 


C m " 4 ! M1 ( p t C t) v 


1 ' = 
C Y = ~ 2tT *M1 (-) Y 


«■ (V) w 


' T 2 \ 1/2 / C 
1 / Z M1 ' 


ir I 


P2 


- Cn \l/2 


( 32 ) 


P ^P 
av o 
n 


Cp i 

$ 


For large magnitudes of control forces and moments for a given power expendi- 
ture, then, as shown in equation (32), the square root of the ratio of ImI 2 to 
Ip2 should be as large as possible. With the general cyclic function of equa- 
tion (27), and equations (28) and (29), values of this ratio for various values 
of m are as follows: 


m 

1 

1/3 

1/5 

0 

3 

2 2 1/2 
( T M1 /! P2 C0S V 

1 

.97 

.95 

.90 

.95 


As can be seen, the sinusoidal cyclic function, m = 1, gives the best value, 
however, this square root ratio is not very sensitive with cyclic function 
since only a ten percent loss is shown for the extreme conditions of a step- 
function, m = 0, cyclic variation. 

These m values, other than one, are rather strong distortions and lead to 
high acceleration points during the cycle. A small variation from sinusoidal 
that uses the next symmetric harmonic is given by 

f(f - y ) = k cos ( 4* - Y ) + (1 - k) cos 3(4' - 4 ' q ) (33) 

which equals the required unity at 4* - 4 'q = 0. Inserting equation (33) into 
equations (24) for n = 2, and (25) for n = 1, results in 


I P2 = k 2 + {1 - k) 2 
'mi = k cos *0 


then the square root ratio is 
1/2 


/ 1 2 ’ 
1 M1 


COS V, 


I 


P2 / 


[1 + ( 


1 


-) ] 1/2 


(34) 


(35) 
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Equation (35) is maximum with k = 1, which is, as before, the sinusoidal 
variation. However, equation (33) will be examined for minimizing angular 
acceleration peaks. Let t = cot, where to is angular velocity, and t is time. 
Taking the second partial derivative with respect to time of equation (33) 
gives the angular acceleration of this cyclic function as 

O 

0 f(Y - y ) 

2 — ~ = <*> 2 |k cos (t - 4* 0 ) + 9(1 - k) cos 3{v - T 0 )J (36) 


The conditions for minima are obtained by differentiating once again and set- 
ting the results equal to zero. Thus 

oj 3 [k sin(i < - h'q) + 27(1 - k) sin 3(y - 4^)] 

(37) 

= cj 3 [81 - 80k + 1 08 ( k - 1) sin 2 (4' - ¥ )] sin(f - v ) = 0 


then 


2 

sin (y 


».> 


80k - 81 _ n _ 81 
1 08( k- 1 ) " u ’ K " 80 


(38) 


The angular acceleration per angular velocity squared for the sinusoidal cycle 
is given by k = 1 in equation (36). With the k of equation (38), the maximum 
angular acceleration per angular velocity squared equals 90 percent of the 
sinusoidal maximum, as shown in the following comparison: 

1 3 2 f(T - Vq) 


T - y 0 , deg 

0 

15 

30 

45 

60 

75 

90 

sinusoidal 

1.000 

.966 

.866 

.707 

.500 

.259 

0 

k - 81/80 

.900 



.899 

.877 

.795 

.619 

.342 

0 


Since vibration fatigue depends on the maximum stress, or acceleration, and 
little on the time rate of change of acceleration, then this k = 81/80 function 
of f^ - M'o)* which is little different from sinusoidal, may be optimum. From 
equation (37) , at v - *Fo = 90 degrees, the k = 81/80 cyclic function has a 
time rate of change of angular acceleration which is 35 percent larger than 
that of the sinusoidal cyclic function. From equations (33) through (35), the 
k = 81/80 cyclic function integral factor characteristics are 
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( 39 ) 


81 1 

“ V = 80 C0S(T " V ' 80 cos 3(f " V 


I = 1 +-§!_ 
X P2 1 3200 


= 1.0253 


81 

J M1 = 80 cos ^o = T- 0125 cos^ 0 
1/2 


1/2 


(I m /I p2 cosS 0 ) = (1 - 1/6562) S 1.000 


The blade angle cyclic variation, cyclic forces and moments, and thrust and 
power, for the k = 81/80 cyclic function, are obtained by inserting equation 
(39) into equations (11) and (18) through (23). Comparing the k = 81/80 
cyclic function characteristics with those due to a sinusoidal cycle (eq. 31), 
the k = 81/80 cycle has the same moment and force ratio per extra power, 1.25 
percent greater moments and forces, 2.53 percent more extra power required, 
and ten percent less maximum blade angular acceleration which helps delay 
fatigue. 

Denotation of Ct linearity with J by CT* . - In equations (18) through (23) 
the partial derivatives of thrust and power with respect to blade angle are 
obtained from steady-state experimental propeller charts or from steady-state 
propeller or rotor theory. Steady-state thrust coefficients tend to become 
nonlinear with advance ratio as advance ratio becomes small. This nonlinear- 
ity is primarily due to the increase in angle of attack on the blade with less 
increase of lift, as that of a wing approaching stall. In addition, the induced 
inflow velocity becomes proportionally larger at smaller advance ratios. The 
blade with cyclic-blade angle sheds plus and minus starting vortices which 
induce velocities that keep the flow on the blade attached, and the added then 
lessened induced inflow velocity during a cycle tends to nullify. These non- 
steady effects thus approximate a thrust linearity with advance ratio as shown 
below. 



With flow remaining attached under cyclic control conditions of no net 
inflow, there would be no drag rise which leads to power increase. Hence the 
partial derivatives with respect to blade angle involving power coefficients 
are taken directly from steady-state experimental data or theory. Taking into 
account the non-steady effect on thrust coefficient, all the derivatives involv- 
ing Ct are with Ct linear with J, denoted by Ct*, for either experimental or 
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theoretical steady-state propeller or rotor data. Then equations (18) through 
(23) become 


1 *" 2 1 °° *Pn 8 n 

Ct - C T = i- I„, C T y 2 + Jr £ ^ y" 


'T ~ °T 4 P2 T 
av o 


. n! 96 " 
n = 4 p o 

even 


1 " 2 1 °° *Pn 9 ^P n 

Co - c n = 1 1 do c D Y + t z Jr- — - - 


' P av ~ " P o " 4 ‘ P2 ^ 


2 n : 4 n! 38 n 
even 


r - r rot u; = It { *r*\' I r n 

C m " C n cot Y o “4 ! M1 ^ p t C T^ Y ' 4 n E ,, n! n Y 

n = 3 9B„ 

odd 


r r rnt v = 1 T f Cp V 1 V If ' 1n n 

C Y -C N cot t q - 2ir I M1 (p p ) Y " 2-rr n f 3 n! ag n Y 

odd 0 


(40) 


(41) 


(42) 


(43) 


where the I integrals are given in equations (24) and (25), and p| is the 
radial position of the center of 1 inear-with-J thrust. 


Cyclic-Control Solutions in Terms of Bessel Functions 

The integrals of equations (2) through (7) for the sinusoidal cycle can be 
directly evaluated in terms of Bessel functions since they already contain 
the integral representations of the Bessel functions. These integrals are 
characterized by trigonometric functions as arguments of trigonometric functions. 
For example, from reference 6 

V 

J (y) = v- sin 2v T cos(y cost) dT, v > - \ (44) 

v r(v + l)r(l) J q 2 

Thus the use of Bessel functions provides a closed solution which is valid to 
large values of the cyclic angle, y. 

Averaged cyclic thrust and power . - The linear with J steady-state thrust 
coefficient is expressed as a Fourier series in b q , thus 


The sinusoidal cyclic-blade angle is given by 


$(* - t q ) = e 0 - y cos(f - t 0 ) 


Substituting b(t - t ) of equation (46) for e Q in equation (45), results in 


c *( f - f Q ) = z s( a n cos ng 0 + b n sin ng 0 )cos[n Y cos(t - t q )] - 
n = 0 L 

(-a n sin ng Q + b n cos ng 0 )sin[ny cos('F - f Q )]^| (47) 

A simple means for obtaining the relationship between integrals of trigonometric 
products and Bessel functions is by use of the following series: 

cos(y sinT) = J q (y) + 2 J 2 (y)cos 2 t + 2J 4 (y)cos 4'? + . . . 

sin(y sin'i') = 2J,(y)sin t + 2J^(y)sin 3 t +2Jr(y)sin 5 t + . . . 

1 J b l (48) 

COS (y COST) = J 0 (y) - 2J 2 (y)C0S 2t + 2J 4 (y)cOS 4 t - . . . 

sin(y cost) = 2J-!(y)cos t - 2J 3 (y)cos 3 t +2J 5 (y)cos 5t - . . . 

and noting that 

•2ir tZ TT -v 

sin mT sin nT dT = cos mT cos nT dT = 0, m f n 


? 2 

sin nT dT = cos nT dT = it, n ^ 0 


sin mT cos nT dT = 0, all m and n 


where m and n are integers. 


The integration of equation (2) with equations (47), (48), and (49), results in 
the average thrust coefficient, thus 
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C* = E (a cos ne + b sin ng )J (n Y ) 
av n = 0 n 0 n 00 


° Ci - E (a cos ng + b sin ng ) [1 - J (n Y )] 
'o n = 1 n 0 n 0 0 

Similarly, for the power coefficient, let 


E (c cos ne 0 + d n sin ne Q ) 
n = 0 


C p = C p - E (c cos ng + d sin ne ) [1 - J (n Y )] (52) 

av K o n = 1 n on o o 

Cyclic pitching and yawing moments . - For integrating equation (4) to 
obtain Cm, the product of the radial center of thrust and the thrust coefficient 
can be represented by another Fourier series. However, since pj does not vary 
much with B 0 it can be approximated as follows: 

p* = pi [1 + k T sin(e - B,)] (53) 

o 


where 


R = _ 3 . /_9B\ 

P 1 4 '9p'p=3/4 


where pf 0 is the radial center of blade, linear-with-J thrust at blade angle 
Bo = 3] , and 3] is the blade angle at p - 3/4 at which the blade pitch ratio, 
p/D = uptan 3 0 (p) is relatively constant with p, that is, pg = constant, then 
B = -p9g/9p. °The constant kj defines the magnitude of the change of radial 
position with blade angle. Substituting Bo(f - To) of equation (46) for Bo 
in equation (53), and factoring this pf function to equation (47), the inte- 
gration of equation (4) results in 


C m = 4 P T COSM, o Z v(- a n sin n 6 0 + b n cos nB Q )J 1 (n Y ) + 
o n = 0 L 

k T (-a n sin n 3 Q + b n cos nB Q )sin(B 0 - B 1 )[J ] (n Y + Y ) + J-|(n Y - Y )]+ 

k T (a n cos nB 0 + b n sin nB Q )cos(B 0 - B-|)[J-|(n Y + Y ) - J-|(n Y - Y )]| (55) 

As an example and a check, let the thrust coefficient be given by Cf = b-| sinB Q , 
then equation (55) reduces to 
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(56) 


c m = j b i pf cosf^EZJ-i (y)cosg n + k T J., (2 Y )sin(2 


r o L " u l 


o T 1 


5 o - P 1 


which is the sinusoidal cycle solution for any value of Y . If the Bessel 
functions are expanded to the third power as 

/M x ) = \ " T6 

then equation (56) becomes 

C m = 1 b l P T C0SY «i c0 SG« + k T sin(2B rt - ejfr - 


o T' 


o P 1 


g- Ecosb o + 4k T sin(2B Q - e-jJJy' 


The check of this example solution is that the power series method of equation 
(42), through n = 3, leads to an identical equation. That is, by using 
f(f - Yo) = ~ cos(4< - 4' ), from equation (46), then I Ml = -1, Im 3 = -3/4, and 
pfcf is the product of equation (53) and b-j sinB Q . 

The integration of equation (5) for C n , is the same as for C m except 
a factor sin4' 0 replaces cos^. Then as before 

C = C tanT n (5) 

n m o 

Cyclic side and normal forces . - For integrating equation (6) to obtain 
Cy, an approximate expression for the radial center of torque or power, Pp, 
can be made as developed for Pf in equation (53). Let 


El + k p sin(3 0 - 


where Pp 0 is the radial center of power at the same conditions cited under 
equation (53), and the constant kp defines the magnitude of the change of 
radial position with blade angle. With the sinusoidal cyclic-blade angle of 
equation (46) inserted into equation (51), the cyclic power variation is 


C (¥ - v ) = i <(c n cos nB 0 + d n sin ng 0 )cosEny cos(r - 4^)] + 
p n = 0 L 


(c n sin nB Q - d n cos nB 0 )sin[n Y cos(¥ - 4^)] 
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Substituting B 0 (t - T 0 ) of equation (46) for Bo in equation (57), and factoring 
this 1/Pp function to equation (58), the integration of equation (6) gives 


COST 

C Y 2irp. 


o 


£ 

n = 0 


|2(-c n sin nB Q + d n cos nB^J^ny) + 


k p(- c n s in nB 0 + d n cos nB 0 )sin(B 0 - B^CJ^ny + y) + J-,(ny - y)]+ 
kp(c n cos nB Q + d n sin nB Q )cos(B 0 - 6-j ) [J-j (ny + y) - J-|(ny - y)]j> (59) 


which, except for notation changes, is the same equation as that for C m in 
equation (55). 

The integration of equation (7) for C^, is the same as for Cy except a 
factor sinT 0 replaces cosTq. Then as before 


C N = -C y tanT 0 (7) 

For the same reasons as those given in the section, Denotation of Cj 
linearity with J by Cf, the steady-state thrust coefficient in equation (45) 
is the Cj that varies linearly with advance ratio J. However, the steady-state 
power coefficient Cp in equation (51) can have a nonlinear variation with 
advance ratio. 


Estimation of Phase Angle 


With cyclic control the blade has a variation of angle of attack during a 
cycle or blade revolution and the blade lift builds up and decreases in a har- 
monic fashion and appears as unsteady motion. As angle of attack increases, 
circulation increases, and a starting vortex is shed which induces a downwash 
and changes the build-up of circulation. In general, unsteady motion causes 
the lift cycle to be out of phase with the angle of attack cycle. An approxi- 
mation for the propeller phase angle can be obtained from an analogous wing 
unsteady solution. In reference 4 an equation is presented which gives the 
phase angle lag of forces and moments due to harmonic pitching of a wing about 
the quarter-chord line, given by 


T 


ph 


k 


= 0.825 tan 

„ %itch 

2V„„ 


2.027k 


(60) 

(61) 


where the phase angle is almost independent of aspect ratio. In equation (61) 
k is the reduced frequency, b is the section chord, and w pitch is pitching 
angular velocity. For the cyclic-control blade the pitching is at the blade 
rotation, then 
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= 2-rrn 


(62) 


“pitch 

The effective velocity the blade senses is Vr as shown below, in which slip- 
stream rotation is taken into account. 



Then for the single propeller or rotor 


1/2 

V R = [(2mr - vtan*) 2 + (V + v) 2 ] , tan* = 2 ^H" 7 ui i ? 

Solving the second equation for vtan<{>, results in 

2 1/2 

vtan<}> = irnr - [(imr) - (V + v)v] 


then with J = V/nD, p =r/R, V R in nondimensional form becomes 


R _ j 1 2 2, , v \ . rl 

HD ~]2 * p + J(J + Hd } + wpc T 17 


2.2 f 1 + v \ v q 1 / 2*1 
nD ; nD J ? 


P c - (J 


1/2 


J 


From momentum theory 


(63) 


(64) 


v 

nD 


= ^[( 0 2 + 




J] 


(65) 


where t is the tip-loss correction, for which approximations developed by L. 
Prandtl are given in references 7 and 8. A modification of Prandtl's correction 
is 


x = 1 


sin 4> 0 
B 




( 66 ) 


which can be evaluated by iteration with equations (65) and (64) for p = 3/4. 
Equation (65) inserted into (64), with p = 3/4, for a single rotor, is 


23 


f 128C. 8 C t ,1/2 

^ 4 + i('-^) ,/2(5 i^ ( ^4) ,/z u (67) 

3 ’ " D /single l. 3 9.V 9/ , T Z J 

For counter-rotating rotors, which are near each other along the axis, w is 
approximately zero. Letting the vtan^ terms be zero in equation (63), then 


. , V + v 

tan * = 2^r 

and Vn with p = 3/4, using equation (65), is 


3tt nD 


counter 


= <1 + 


32C t 8C t 

o + 8 w t (J 2 + _^ ) , /2 ] 


where Cj Q is the thrust coefficient of one of the two counter-rotating rotors, 
or one-half of the CT 0 of the total propeller. From equations (61) and (62), 
the reduced frequency is 


4 3/4 

3 D 

4_ \ 
3ir nD 


3/4 o 


4 1R 
3tt nD 


where B is the number of blades on each rotor, and a is the effective solidity 
of each rotor, defined in the relation 

a=^-V (71) 

3 TT D 

The phase angle from equation (60), with equation (70), is 


6.368 


f. = 0.825 tan" 
ph 


3/4 a 
b av 6 


4_!a 

3tt nD 


where the denominator in the arc tangent function is given in equation (67) for 
the single rotor, and given in equation (69) for each of the counterrotating 
rotors. The phase angle is the azimuth angle lag which thrust and torque lag 
the blade angle variation. Thus, maximum thrust occurs an angle after maxi- 
mum blade angle. The basic azimuth angle, which governs the components of 
the cyclic-control forces and moments, is given in equation (1). By using 
equation (1), for counter-rotating rotors, the clockwise rotation rotor has. 
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1 


Y q = ^ph + ^ax + ^con’ w * 1 ^ e the counterclockwise rotation rotor has, 

Y 0 = “^ph + + ^cop- The an 9^ e can be ma( ^ e the same f° r the counter- 

rotating propeller by using a different control angle, Ycon> on each rotor. 
Since the variation of V^/nD with thrust coefficient is small, the phase 
angle will change little with thrust. The Vr velocity increases moderately 
with advance ratio, then the phase angle will decrease as J increases. At 
static conditions, J = 0, the phase angle will slightly increase with thrust 
increase for the single rotor, and decrease with thrust increase for each of 
the counterrotating rotors. 

This derivation of phase angle is based on the wing, or blade, with non- 
separated flow. If the cyclic-blade pitching is made about the stall C£ of 
the blade, then it is operating in a C£ ~ a hysteresis loop in which the air 
adds energy to the blade. Then the phase angle is approximately given by a 
negative value of equation (72) when k < 0.3, as shown by an experimental 
example in reference 9. 


CYCLIC-CONTROL COMPARISONS OF THEORY WITH EXPERIMENT 
AND EVALUATION OF EXTRA POWER 


The cyclic-control theoretical results of equations (40) through (43), 
and (50), (52), (55), and (59), are based on knowing the steady-state values 
of thrust, thrust linear with advance ratio, and power coefficients. Here 
steady state refers to the condition that maximum cyclic-blade angle is zero, 
y = 0, and Cj, Cj, and Cp versus J for the given propeller are obtained from 
available propeller chart data, or from a favored propeller theory. An approx- 
imate steady-state prediction method from reference 4, for estimating thrust 
and power coefficients and their derivatives, is further developed and pre- 
sented in appendix B. The experimental study reported in reference 1 contains 
propeller cycl ic-control experimental data and steady-state data for J = 0. 

The objectives are to show an example usage of the theory, correlate with 
experiment, and develop relations for extra power due to cyclic control. The 
cyclic-control power is the ideal power plus extra power where extra is that 
power other than ideal traceable to the varying cyclic-blade angle. 


Comparison of Theory with Experiment 

When the cyclic-blade angle y is small, then the terms y 11 for n > 2 
are negligible. Then equations (40) through (43) for the sinusoidal cycle, 
f ( Y - Yo) = cos(y - Yq), simplify to 


C 






(73) 
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'Cp -Cr 
av 

2 r 

, Y Cp 


1 s_ 

4 Cr 


( 74 ) 


ideal 


m 


C* C* 


cr ' C 1 cot ’o = - 4 (p T t~ + 


)cos V 


P T ‘0 


(75) 


J N 


p— ^ = - -p-X- cot ¥=-*-[— p“— + (—) ]cos v 

Cn Cn 0 2 tt pp Cp Pp 0 


i 


(76) 


In the experiment, all the quantities on the left side of equations (73) 
through (76) are measured at the condition of J = 0. Also the steady state, 
that is, Y = 0, thrust and power coefficient are measured at J = 0. In this 
section the experimental steady-state values of Cj and Cp are used to determine 
the thrust and power derivatives. Using these derivatives, the theoretical 
values of cyclic-control forces, moments, change in thrust, and change in power, 
are determined from equations (73) through (76), and compared with measured 
values. 

Cyclic-control moment and change in thrust . - The experiment values of Ct 
were measured at J = 0. For the cyclic-control theory a relation between Cj 
and Cf is that in equation (B18) of appendix B. Since both 9CT/8J and JoT in 
equation (B18) are functions of 3 p , then the Cf derivatives are obtainable. An 
approximate Cj ~ Cf relation is given in equation (BIO), and Cf derivatives in 
terms of Cj derivatives are given in equations (B19) through (B22). These 
derivatives are applicable for the J = 0, Cy measurements. 

The geometric characteristics of the propeller used in the test are as 
follows : 


B = 3, — jp = .1178, a = .15, = 3 + e, e = 2° 

J = 0, k = .9, <5g = .015, $1 = 17°, d Q = .025 ^ (77) 

e T = .343, e p = .634, J qT = 2.246 tan 3 Q 

where a, 3-| , d 0 > ej, e p , and J 0 j are determined from equations (71), (54), (B9), 
(B3), (B4) , and (B5), respectively. Curve fitting a function for the experimental 
static steady-state thrust data, leads to the relationship 

C T = .296 sin (2 b + 3°) (78) 
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then the derivatives with respect to 3 are 


T = 2 

Cj _ tan(2g + 3 C 


= -4 


An approximation of the pf function is given in equation (53). 
vative with respect to 3 0 is 

p*' = P f 0 k T cos(e 0 - 6,) 

For typical propellers 


The deri- 


p To - 3 ’ k T =7 sin + 5°) 1 

With the values of equation (77), and with equations (53), (81), and (82), 
equations (73) and (75) become 

r - r » 

aC T T av T 1 C T 
2 r 2 r 4 C T 


C T cos ¥ 0 = ~ 6 i [1 + - 381 Sin (e o 


cos(3 q - 17°) 


17°)] + .381 E 1/2 


★ 1 I I 

where Ct /Cj, and Cj /Cj, from equations (B19) and (B20), with equations (79) 
and (80), are 


+ n_] E V 2 

E T J h T 
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if 


i 


t 


ii 


E T 1 , E T \Z . 1 e t 


1/2 


T _ r A + 2 T I / T J_ , 1 F 

“CjT tan (2g + 3) E-j. " 4 ' 2 E ? JE T 

where from equations (B12), (B21), and (B22), Ej and derivatives are 


( 86 ) 


Ej = 1 + .7353 (1 + .375 sin 26 Q )‘ 


Ej ° 1.1030 (cos 2 b q + .1875 sin 4 b q ) 


Ey = 2.2059 (.375 cos 4 b q - sin 23 Q ) 


(87) 

( 88 ) 
(89) 


The values determined from equations (78), (83), and (84), are compared 
with experiment for given values of blade angle, as follows: 


e 

deg 

e o 

deg 

C T 

AC t /y 2 C t 

-Cm Y /C T cos 

eq. 

(78) 

experi- 

ment 

eq. 

(83) 

eq. 

(84) 

experi - 
ment 

8 

10 

.0963 

.092 

- .213 

1.438 

1.33 

10 

12 

.1156 

.112 

- .493 

1.222 

1.22 

12 

14 

.1344 

.134 

- .708 

1.066 

1.03 

14 

16 

.1525 

.153 

- .808 

.946 

.93 

16 

18 

.1698 

.172 

-1.024 

.849 

.85 

18 

20 

.1863 

.188 

-1 . 1 68 

.767 

.78 

20 

22 

.2019 

.203 

-1.256 

.696 

.72 

22 

24 

.2164 

.216 

-1.350 

.632 

.66 


As shown in these results, the Cj of equation (78) duplicates Cj of the experi- 
ment quite accurately. For y not large, the theoretical thrust loss is negli- 
gible at small 3 , but with greater loss at large blade angles. In the experi- 
ment the thrust has a slight loss only at the higher 3 angle. The pitching 
moment coefficients correlate within the accuracy of the measurements. 

Cyclic-control force and ideal change in power . - Development of an analy- 
tical function for a curve fitting the experimental static steady-state propeller 
power data, leads to the relationship ; 
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(90) 


C = .292 tan (8 + 11P ) 

P cos 3 B 

then the derivatives with respect to $ are 


CT = tan ~ l ' g 2 ' + T P ' ) ' + 2ta " < B + n °> + 3 tanB 


(91) 


c c 

-£ — = 5 — + 11 + 6 tan 3 (g + 11°) - 6 tan^g + 6 tanB (92) 

L P tan^ (b + 11°) P 

where b = B 0 _ 2 °. 

The function 1/Pp is needed in equation (76). An approximation of 1/Pp 
is given in equation (57). The derivative with respect to Bo is 


( 7 -) = ~~ cos (p - Bi ) 

Pp Pp 0 1 


(93) 


For typical propellers 


Pd - n > k n £ - 


Q " 11 * K P = " 11 sin (b ] + 5°) 


(94) 


With the propeller values listed in equation (77), and with equations (57), (93), 
and (94), equations (74) and (76) become 


AC r 


ideal _ 1 *av 


2 r 
Y C p 


y C c 


1 Si 

4 C n 


(95) 


/ideal 


C p cos 


= -.2188 <[1 - .243 sin (r - 17°)] - .243 cos 


(s 0 - ,7 °>) 


(96) 


where Cp, C p ', and Cp ' 1 are given in equations (90), (91), and (92), respectively. 
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The values determined from equations (90), (95), and (96) are compared 
with experiment for given values of blade angle, as follows: 


6 

deg 

6 o 

deg 

c p 

p /T 2 C p 

-CY/C 

o 

l/> 

-6 

O 

eq. 

(90) 

experi- 

ment 

eq. 

(95) 

experi - 
ment 

eq. 

(96) 

experi- 

ment 

8 

10 

.0357 

.036 

8.574 

16.1 

1.505 

1.17 

10 

12 

.0451 

.046 

8.039 

13.3 

1.400 

1.28 

12 

14 

.0563 

.057 

7.704 

11.6 

1.320 

1.35 

14 

16 

.0695 

.070 

7.515 

9.8 

1.259 

1.44 

16 

18 

.0853 

.084 

7.439 

10.7 

1.211 

1.15 

18 

20 

.1043 

.102 

7.454 

9.7 

1.176 

1.00 

20 

22 

.1270 

.127 

7.547 

10.5 

1.151 

1.21 

22 

24 

.1546 

.155 

7.711 

8.7 

1.134 

1 .00 


It can be seen that the Cp of equation (90) duplicates Cp of the experiment 
quite accurately. The experimental data of side force coefficient is somewhat 
erratic, however, the predictions show a correlation with experiment of magni- 
tude and that CYy/Cp decreases with increasing e. The predicted change in power 
due to cyclic control is the ideal change in power, analogous to induced drag. 

The obviously higher measured change in power includes the ideal power change, 
plus extra or additional power expenditures. This extra power includes power 
required effects peculiar to cycle control but nonexistent in steady-state 
operation. Extra power characteristics are investigated in a subsequent section. 


Extra Power Changes Due to Cyclic Control 

Extra power changes include energies due to a vibratory swirl of the 
slipstream, angled slipstream, pitching airfoil of the propeller blades, and 
vibratory structural damping. An estimate of these extra power expenditures 
can be made without excessive complication. 

Vibratory swirl of the slipstream . - For the steady-state propeller, part 
of the power is used to swirl or rotate the slipstream, and this part equals the 
time rate of change of the kinetic energy of the slipstream rotation. Similarly, 
the cyclic-control propeller swirls the slipstream in one direction, then rever- 
ses the swirl during the other half of the cycle, or results in a vibratory 
frequency equal to the propeller revolutions per second. Thus the cyclic-control 
propeller leaves vibratory swirl kinetic energy in the slipstream which the 
steady-state propeller does not. This kinetic energy per unit length of the 
slipstream is 

K.E . = 1 1 • 2 

£ 2 £ sw 
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where the moment of inertia per unit length of a circular column of air is given 
by 

? 

I mass R _ it n 4 
z z 2 2 p a K 

and is the swirl angular velocity. With cyclic control, Y S w = ^swo cos tat, 
then v sw = -'i’swo a) sin cat. The kinetic energy in a complete revolution is the 
integration of *2 SW> or 


K.E. _ 1 I 


where Tswo is the maximum swirl angle due to y. The swirl power is the kinetic 
energy per unit length times the velocity increase through the propeller disc, 
then with the I/z expression and since oj = 2trn 


4o " 2 s1n2fd * 


It I 2 2 

21“ *swo 


P 


sw 


4 2 4 2 
p a ir n RV 


swo 


V 


or 



= ll ^ v_ 

16 swo nD 


(97) 


where v is the average velocity increase through the propeller disc. From 
simple momentum theory, for example, reference 7, p. 183, the inflow velocity 
in terms of Cp is obtained by the solution of a cubic equation, then 


v _ 2 rJ 3 C P ( 2J % C p 2 ) l/2 1 l/3 rJ 3 C P 

nD 3 J + [ 27 + r + { ~ThT + ^2> ] + [ 27 + ~ 


•27 


27tt 


- (- 


2J 3 C p 
27F 


^_) 1 / 2 ] 1/ 3 


J = 0 2C, 


- (- 


P)l/3 


(98) 


Similarly, the maximum swirl angle can be estimated from momentum theory. Then 


<' ’ 5"'swo )( ? ’swo' 


2C t p = 2/3 9 C t 


3 2 

7T p 


2tt' 


where here Cj = cf(So + y) - Cf(g 0 )- This Cj can be approximated by using equa- 
tion (Bl), then 

9e 

T swo S ~r (2 - 2 coss o + J sine o )Y (99) 
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With equation (99), the power expended for vibratory swirling the slipstream, 
equation (97), becomes 


if = h ( tt ) 2 {2,2 cos % + J sin6 o )2 HD 


(100) 


Y C 


where v/nD is given in equation (98), and ej in equation (B3). For J - 0, 
equation (100) simplifies to 


P e T o 

SW _ o ioc T 

2„ 2,136 r 2/3 C0S 3 o 


(101) 


Y C, 


Angled slipstream . - When the cyclic-blade angle is applied for pitching, a 
side force is developed. The reaction to this side force produces a slipstream 
angle opposing the side force. This slipstream angle, represented by e Y » equals 


_1 c y C Y 

e Y = tan = r~ 

y c T c T 


or for an arbitrary angle of f 0 > and using equation (7) 

C v 


e = 
Y 


- ^ (Cy cosf 0 - C N sinT 0 ) cocos' 


( 102 ) 


Because of the conservation of momentum, due to slipstream constriction the 
slipstream velocity increases approximately as 

V- 

V . ( e ) = , where = V + v 
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From momentum theory 
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The ratio gives 
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with equation (102) 
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With these relations, the change in power coefficient, C p - Cp , is the angled 
slipstream power coefficient, given by e 




C T cosy 
T o 


(1 



(103) 


where v/nD is given in equation (98). Equation (103) is an estimate of the 
fractional power increase due to the angling of the slipstream due to cyclic 
angle y. It can be noted that for the counterrotating propeller if one rotor 
operates at y 0 , and the other rotor at it - then by equation (82) the net 
Cy is zero. Thus from equation (103), Cp is zero, and from equation (102), 

e is zero. y 

Y 


Pitching airfoil of the propeller blades . - Energy is imparted to the 
airstream due to forced pitching of an airfoil. An estimate of this energy per 
unit time is given in reference 9, page 169. For pitching airfoil blades, the 
power due to blade pitching is 


IT 

16 



3 2 

w y 


where in coefficient form is 



Y 


4 


where b gv /D is related with effective solidity in equation (71 ) . 


(104) 
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Vibratory structural damping . - The power absorbed in structural damping 
in a propeller is small, but can be larger for the less rigid helicopter blade. 
Structural damping is discussed in reference 9, page 227. The power expendi- 
ture for damping, due to pitching, is approximated as 


r2ir 


damp 


BI g a) 

a ot ot 
2io 


V 3t ; 


dv = 


77 RT n 2 2 

7T BI g a) wy 
2 ora a 


where I a is the moment of inertia of the blade, is the natural angular 
frequency of the blade, and g a is the damping coefficient which is a fraction 
of the elastic coefficients. Empirically, g a is found to depend on the ampli- 
tude of motion, and power could probably be better represented as proportional 
to y2. The moment of inertia of a blade is approximately 


I = 


a 


TT 
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P 

b 


t(l + t 2 )b 4 D 


where pb is the blade density, t is maximum thickness to chord ratio, b is 
blade chord, and D is propeller diameter. Then the structural damping power 
coefficient is 


- P ltm , l! B(-jp) 4 t(l + t 2 )(V ^ g 

32 D a) p_ a 


y 2 ~ "a 

where Pb/pa is the ratio of blade density to air density. 
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Total Power Comparisons of Theory with Experiment 


The total power is the sum of the ideal and the four extra power require- 
ments of the cycl ic-control propeller or rotor. The test propeller character- 
istics are listed in equation (77). In addition, estimates of the blade 
density ratio, natural angular frequency of the blade, blade damping coefficient, 
and thickness ratio are 


(— ) 2 = 4000 

p b 

— = 2200 
p a 

g = io " 5 

a 

t = .09 


(106) 
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which approximate that of an aluminum blade with tip speed of about .8 Mach 
number. The damping coefficient of a wing is typically of the order g a = 10-3 
to 10 - 2 . The propeller blade is much more rigid than a wing and g a is assumed 
to be one percent that of a wing. For these propeller values, the power coef- 
ficient changes due to cyclic-control are compared in the following table: 


ACp/y^Cp 


0 

deg 

ideal 

eq. 

(95) 

swi rl 
eq. 
(101) 

angled 

eq. 

(103) 

pitch. 

eq. 

(104) 

struc. 

eq. 

(105) 

theory 
total , 
5-part 

experiment 

8 

8.57 

2.61 

.44 

.40 

1 .20 

13.22 

16.1 

10 

8.04 

2.23 

.42 

.31 

.93 

11.93 

13.3 

12 

7.70 

1.92 

.42 

.25 

.75 

11.04 

11.6 

14 

7.52 

1.67 

.45 

.20 

.60 

10.44 

9.8 

16 

7.44 

1.46 

.50 

.16 

.48 

10.04 

10.7 

18 

7.45 

1.27 

.58 

.13 

.39 

9.82 

9.7 

20 

7.55 

1.12 

.68 

.11 

.33 

9.79 

10.5 

22 

7.71 

.97 

.86 

.09 

_ 

.27 

9.90 

8.7 


The proportion of power expenditure due to vibratory-swirl is relatively 
large and it appears this could be reduced by synchronized counterrotating 
blades, that is, blades positioned so that maximum swirl is nullified. 


COUNTERROTATING PROPELLERS AND OTHER CONFIGURATIONS FOR 
ISOLATING CYCLIC-CONTROL FORCES AND MOMENTS 


In any control system it is desirable to initiate a given force or moment 
without other accompanying forces and moments. For example, to initiate pitch- 
ing moment without side force or thrust change, that is, to isolate pitching 
moment from the other control forces and moments. 

A single propeller with cycl ic- control generates 
a moment and a force and neither can be isolated. 

As shown here schematically, a positive y applied 
to gain pitching, also gains a side force. The 
side force comes from the reaction to the 
increased torque on the blade when the blade 
is at a e 0 + y angle at ¥ = it. The + sign 
indicates a positive moment where y maximizes 
the blade angle. This single propeller effect, of forces and moments which 
cannot be isolated, can be seen analytically in equations (75) and (76). For 
any value of azimuth angle in these equations, there will be a force and moment 


propeller rotation 

v^-aft view of 
rf — ^ \ propeller 
side V i disc 
force y&S 

M, moment, y at T = n 
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combination. However, the counterrotating propeller has unique cyclic-control 
qualities of isolating control forces and moments, and in addition, rolling 
moment control, that is, moments about the axis. Multi-propel lers , in-line, 
can also isolate forces and moments. 


Counterrotating Propeller 

Two blades of the counterrotating propeller with cyclic-control, at two 
different operating modes, are shown in aft view as follows: 



In the left drawing, with a plus cyclic-angle yon one blade and a negative 
angle yon the counterrotating blade, the moments cancel, while the reaction 
forces add and provide a net translation force. In the right drawing, a plus 
angle yon both the downgoing blade and the counterrotating blade results in a 
cancellation of forces and the moments add to provide a net moment on the 
propeller plane. Since azimuth angle Y is not specified in this example, then 
the results apply at any azimuth angle, hence with an azimuth control angle, 
forces and moments can be controlled singularly or jointly in any direction 
in the propeller plane. Thus the cyclic-control counterrotating propeller has 
a built-in capability of isolating all forces and moments or combinations, 
without the need of propeller tilt or of net thrust changes. The above left 
drawing can be considered as a rotor in a horizontal plane, that is, helicopter 
mode, then the cyclic-control force becomes a direction-controllable propulsive 
force while the rotor remains horizontal, or untilted. For this configuration 
the rolling moment, that is, a moment about the propeller axis or shaft, is 
equal to a torque differential made by a blade angle difference between the 
two rotors. 

Cyclic-control analysis . - The analysis developed for the single propeller 
can be applied to the counterrotating propeller by taking into account the azi- 
muth angle for control. From equation (1) it can be seen that the azimuth con- 
trol angle is linearly proportional to Yo, hence, y 0 will be used as a displaced 
control angle. The cyclic-control solutions with sinusoidal cycle, f(Y - Yo) 

= cos(y - y 0 ), for the single propeller are given in equations (73) through 
(76). For the counterrotating propeller let subscript 1 refer to the rotor 
with clockwise rotation, and subscript 2 to the rotor with counterclockwise 
rotation. Then equation (73) applied to each of the two rotors becomes 
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Adding these two equations gives 
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For equal and equal operating clockwise and counterclockwise rotating rotors 


* 1 1 * 1 1 * i > 

C T1 = C T2 = C T /2 


Then 
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(108) 


where these thrust coefficients and derivatives are of the counterrotating pro- 
peller. Similarly, for the power change 


(On - Cp) « ui '-pi 

K av v ideal 4 1 
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When Y i = Y 2 = Y > equations (108) and (110) become the same equations as those 
for the single propeller. 


From equations (75) and (76), define K’s as 

p * i p * 

v _ 1 ( * T , U T * n 

\\ 4 P T C T C T P T ^ 


(111) 




(112) 


From equation (75), for each of the rotors 

C ml = _C T1 K ml Y 1 cos ' i '01 ’ C m2 = " C T2 K m2 Y 2 cos ' 1 '02 
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Then adding 


C m "£ C Tl K ml Y l COS ' 1 'oi 

+ C T2 K m2 Y 2 COS4 '02-' 

(113) 

Similarly, then 



C n = -[ C Tl K ral''l sinT 01 

+ C T2 K m2 Y 2 sinvF 02^ 

(114) 

With the same procedure, using equation (76) 


C Y “ t c pi k yi ^l cos ' i 'oi 

- Cp 2 Ky 2 Y2 ( -02 V Fq 2 D 

(115) 

C N = C Pl K Yl Y l siny 01 

C P2 K Y2 Y 2 Sin4 '02 

(116) 

where the negative signs of the second terms are because the torque reaction 
force of counterclockwise rotation is in opposite direction to that with clock- 
wise rotation. For equal and equal operating clockwise and counterclockwise 
rotation rotors 

Cyi Cy2 Cy/2 , 

K ml = K m2 = K m’ 

C P1 C P2 C P /2 ] 

K Y1 = K Y2 K Y J 

> (117) 

then equations (113) through (116) become 
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For isolating the forces from the moments, the terms within the parenthesis | 
of equations (118) and (119), must be zero. Then 

I 

y 2 -COS¥ 01 _ COs(tt + H'qi ) Y 2 _sin ^01 sin ^ + W 
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Equating these two equations, ¥02 can be evaluated as 


^02 “ 71 + T 01 ’ then Y 2 = Y 1 022 ) 

Similarly, for isolating the moments from the forces, the values within the 
parentheses of equations (120) and (121) must be zero. The solution gives 

T q2 = ¥qi. then y 2 = Y] (123) 

Thus, when the cyclic-blade angle magnitude is the same on both rotors, the 
cyclic-control forces and moments can be isolated with a counterrotating pro- 
peller by the azimuth control angle on each rotor. With y-j = Y 2 = Y, equations 
(118) through (121) become 



= 4 V cos, 01 + cos ' t 02 ) 
= 4 K m (si, "01 * s1nl * , 02 ) 

= Ky ( cos^qi - cosv 02 ) 
= \ Ky(sin^ 01 - sin¥ 02 ) 


024) 


025) 
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where Tqi is the To of the rotor with clockwise rotation, and T02 is the To of 
the rotor with counterclockwise rotation. The K's are given in equations (111) 
and (112). 


The rolling moment coefficient about the propeller axis or shaft is 

C r = — 4 t 5 = 4r [C P (6 o + Y l> ' C P (e o ' Y 2 )] (128) 

p a n 0 


where the blade angle of rotor 1 is increased by the angle y-j, and of rotor 2 
decreased by Y 2 , where L is the rolling moment, and the parentheses represent 
functional parameters of the power coefficient. Equation (128) is the change in 
torque between one rotor and the other, then 
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= -k [C P n (3 o + Y l } ' C P„ (3 o ~ Y 2 )] 


which with Cp] = CP 2 = Cp/2 for equal rotors, leads to equation (128). 


Equations (124) through (127) are such that any arbitrary control force, 
moment, or combination can be initiated by the azimuth angle control of *oi and 
^02 • For the case of isolated moments, Tq 2 = ^Ol > then 
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The resultant moment is 


C (C 2 + C 2 )^ 2 

mn ' m L n ' 

I - Y Y 


= Km’ at "mn = "oi 


(130) 


where y mn is from tan% n = Cny/Cmy. This resultant moment is acting forward 
and normal to the propeller disc at azimuth angle ^oi ■ For the case of 
isolated forces, ^02 = ir + ^01. then from equations (124) through (127) 
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The resultant force is 
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where ^yn is from tan = C|s| y /C Y . This resultant force is acting radially 
outward from the shaft at azimuth angle '•'yn* Thus to isolate moments, Tq 2 1S 
set equal to ^oi » and to isolate forces, * 1 'o2 is set equal to ir + ^Ol > and 
TqI is controllable through the T C on control angle in equation (1). 


The cyclic-control analysis for the counterrotating propeller include 
equations (108) and (110) of the change in thrust and power, equations (124) 
through (127) of cyclic-control forces and moments, equation (128) of rolling 
moment by blade angle differential, and equations (130) and (132) of isolated 
cyclic-control moment and force. The Km and K Y functions are given in equations 
(111) and (112). The extra power required in cyclic-control operation, approxi- 
mated in equations (100), (103), (104), and (105), apply also to the counter- 
rotating propeller, except for the isolated moment condition. For this condition 
the vibratory swirl from the two rotors cancel each other and equation (100) is 
assumed zero. In the isolated moment condition, C Y = Cn = 0, then equation 
(103), angled slipstream power, becomes zero. Since force along the propeller 
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axis, that is, thrust, is simply controlled by steady-state blade angle, then 
six degrees of freedom are controllable with a cyclic-control counterrotating 
propeller. 

Unique cyclic-control characteristics . - The counterrotating propeller 
capability for isolating six control forces and moments is particularly ideal 
for large lifting-crane type vehicles, for airship and submarine control, and 
for tilt rotors. The lifting-crane helicopter can be made to generate a pro- 
pulsive force in any horizontal direction by isolating the side and normal 
forces, done without tilting the rotor. For hovering, the six degree of freedom 
direct control can be maintained without any movement of the aircraft, making 
it an exceptionally stable platform. For the same reasons, a cyclic-control 
counterrotating propeller, concentric at the tail of an airship or submarine, 
can control the hover of these vehicles without vehicle movement, or at forward 
velocities can provide six degree of freedom direct control. 


Multipropellers In Line 


Single rotation propellers when in line as on a planar wing, can be 
operated so that some or all the control forces and moments can be isolated. 
Examination of control forces and moments can be made from drawings of the 
propeller disc, represented by a circle depicting an aft view of the disc. The 
curved arrow on top of the disc indicates the direction of propeller rotation, 
the straight arrow indicates control force, the + sign shows position of cyclic- 
angle y, and the control moment acts out of the page at azimuth angle position 
denoted by +. Results for two propellers in line are as follows: 

1 2 


to isolate with equal y, side forces 

pitching: are cancelled. 


yawing and 

rolling: cannot be isolated. 



normal 

force: with equal y, yawing and 

rolling are cancelled 



side 

force: with equal y, pitching is 

cancel led 
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Hence, for two propellers, all control forces and moments can be isolated except 
yawing from rolling. 


Results for three propellers in line, are 


to isolate 

pitching: with = y-j + y 3 » y-j = y 3 > 

side force and yawing cancel 


yawing: with y£ = y-j + y 3 » y-j = y 3 

normal force and rolling 
cancel 

rolling: with Y 2 = Y-j = y 3 

normal force and yawing 
cancel 


normal 

force: with Y 3 = y-j + y 3 * y-j = Y 3 

yawing and rolling cancel 


side 

force: with Y 2 = y-j + y 3 

pitching is cancelled 


1 2 3 






With three propellers, propeller rotation given by the sense \ /""'>( 

can also isolate control forces and moments, however, the relationship between, 
y 1> y 2> y 3 for isolating rolling, becomes complicated; also, the yawing ability 
is weak since it would be only that due to the center propeller moment arm off- 
set. 


Results for four propellers in line, are 

to isolate 
pitching: 

yawing : 
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rolling: 


normal 

force: 



side 

force: 



With four in line propellers, the rotation sequence, ^4 ~n , of 
the group on the left is preferred because the isolated yawing moment magnitude 
is twice that obtained from the rotation of the group on the right. The other 
isolated control forces and moments are the same in either group. 


In all of the above configurations, an opposite propeller rotation leads 
to the same result, except the directions of the straight arrows are opposite, 
leading to some sign changes in the forces and moments. It can be noted that 
three is the minimum number of propellers grouped laterally which can, by 
cyclic-control, isolate all the control forces and moments. Cyclic-control 
analysis for each of the propellers in a multipropeller group is assumed to be 
the same as for the single propeller analysis given in equations (73) through 
(76). Greater accuracy is possible if the steady-state Cj and Cp versus blade 
angle values are known for each propeller in a multipropeller group. 


RESULTS AND DISCUSSION 


Guide for Usage of Theory 

The theoretical cyclic-control equations are derived as a generalized theory 
relating cyclic-control forces, moments, power, and thrust changes, to the 
steady-state characteristics of the same propeller. For each cyclic-controlled 
propeller there are two arbitrary controls. These include the cyclic-blade 
angle y, for governing the magnitude of the cyclic-control force or moment, and 
the azimuth control angle, ^on* for governing the direction or sense of the 
cyclic-control force or moment. 

Sinusoidal cycle . - For the cyclic-blade angle which varies sinusoidally 
with azimuth angle, the cyclic function is f('F - f 0 ) = cos(f - ¥ 0 ), then with 
angle y not large, the cyclic-control relations are given in equations (73) 
through (76). These relations are most useful, and are theoretical ly accurate 
but analytically simple. Equations (40) through (43) apply for the condition 
of large y. For the sinusoidal cycle the Ip n and Iw n are evaluated by substi- 
tuting n for a in equations (28) and (29), then iP n = 1, 3/4, 5/8, 35/64, . . ., 
(n - l)Ip,n-2/n, n = 2, 4, 6, . . ., even, and Ipin/cosTo = 1 , 3/4, 5/8, 35/64, 

. . ., n I m 5 n _ 2 / ( n + 1 )cosyq, n = 1, 3, 5, . . ., odd. 
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The cyclic-control solution in terms of Bessel functions represents a 
closed function solution for any value of y. This solution is in terms of the 
Fourier coefficients of the steady-state thrust and power and radial centers 
of thrust and power. The linear-with-J steady-state thrust coefficient is 
represented in a Fourier series in Bo as in equation (45). Then the cyclic- 
control change in thrust is given in equation (50), and the moments Cm and 
Cn are given in equations (55) and (5). Similarly, with the steady-state power 
coefficient represented in equation (51) as a Fourier series, the cyclic-control 
ideal change in power is given in equation (52), and the forces Cy and Cn are 
given in equations (59) and (7). The blade angle variation in a cycle is given 
in equation (46) where Yo from equation (1) includes the azimuth control angle 
fcon- 


An alternative solution in terms of Bessel functions is as follows. If 
the steady-state theoretical or experimental data includes the blade radial 
center of thrust, then the product pfCf can be expressed as a Fourier series, 
given as 

oo 

p*ci = E (e cos ne + f sin n3 n ) (133) 

11 n = 0 

where e n and f n are Fourier coefficients. The cyclic-control solution for C m 
is the same as that of equation (55) except en and fn replace an and bn> and 
since py is part of equation (133), Ky = 0, py = 1» then equation (55) becomes 

oo 

C m = C n cot'i ' 0 = ^ C0SY 0 E (-e n sin nB Q + f n cos nB^J^ny) (134) 

Similarly, for Cp/pp expanded in a Fourier series 

C 00 

7T = E (9 n C0S n ^n + h n Sin nB n ) ( 135 ) 

p P n = 0 

then with equation (59) 

1 

Cy = -C^cotfj = -7 cos ¥ e (-g n sin nB Q + h n cos nB 0 )J-j(ny) (136) 

Steady-state propeller data . - Propeller theoretical or experimental per- 
formance charts contain Cy and Cp as functions of advance ratio. From the 
Cy ' J chart, Cy from equation (B18) can be evaluated as a function of steady- 
state blade angle, B 0 . With the Cf(Bo) function known, the derivatives of 
Cy with respect to 3o> or Cf as a Fourier series in 3q, can be formulated for 
use in the cyclic-control theory. Similarly, for a given value of J, the 
Cp(3 0 ) function is known, then the derivatives of Cp with respect to 3o> or 
Cp as a Fourier series in 3o> can be formulated. An approximation for the 
radial center of thrust is given in equation (53) and its derivative in equa- 
tion (81). A similar approximate expression for radial center of power or 


44 



torque is given in equation (57) and its derivative in equation (93). Propeller 
performance data, particularly theoretical, can include the change of thrust and 
power with change in radius, as a function of blade radius position. Then the 
radial centers of thrust and power as functions of 3 0 given directly for 
application in the cyclic-control solution equations. 

An alternative but approximate steady-state solution is presented in 
appendix B. The first and second partial derivatives of the thrust coefficient 
are given in equations ( Bl 6) and (B17). Equation (B17) into (73) gives the 
cycl ic-control thrust change, and equations (B16), (BIO), (53), with (75) gives 
the cyclic-control moments. The first and second partial derivatives of the 
power coefficient are given in equations (B23) and (B24), which when J = 0, 
simplify to those in equations (B31) and (B32). With equation (57), the ideal 
power change and cyclic-control forces are determined from equations (74) and 
(76). 


Arbitrary cyclic function . - The cyclic function is f(4» - ^ 0 ) defined in 
equation (11). The cyclic function can be any arbitrary antisymmetric function 
with azimuth angle, but by definition, bounded within ±1. The generalized 
cyclic-control solution for an arbitrary cyclic function is given in equations 
(40) through (43). Example solutions of arbitrary cyclic functions are given 
in the section, Cyclic function of cosine to m-power. 

Extra power changes . - The total power change due to cyclic-control is 
the sum of ideal power change plus extra power changes which include vibratory 
swirl of the slipstream, angled slipstream, pitching airfoil of the propeller 
blades, and vibratory structural damping. An estimate of these extra power 
expenditures can be made from equations (100), (103), (104), and (105), respect- 
ively. 

Counterrotating propeller . - With the condition that each rotor has the 
same value of Y, the cyclic-control pitching and yawing moments are given in 
equations (124) and (125) and side and normal forces in equations (126) and 
(127). When moments are isolated, that is, forces are zero, the resultant 
moment is given in equation (130) at a given azimuth angle. When forces are 
isolated, that is, moments are zero, the resultant force is given in equation 
(132) directed radial outward from the shaft at a given azimuth angle. The 
thrust change is given in equation (108), and the ideal power change in equation 
(110). The rolling moment about the propeller shaft is the change of torque 
between the two rotors, given in equation (128). The longitudinal force along 
the propeller axis, or thrust, is simply controlled by steady-state blade angle. 
The extra power required in cyclic-control operation is approximated in equations 
(100), (103), (104), and (105). However, in the isolated moment condition the 
vibratory swirl from the two rotors cancel each other and equation (100) is 
assumed zero. Also in the isolated moment condition, the forces are zero, then 
equation (103) becomes zero. 
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Example Solution and Experimental Comparison 

An example use of cyclic-control theory is shown by applying the theory to 
a test-model propeller. The theory is given in equations (73) through (76). 

The geometric characteristics of the propeller are listed in equation (77). 
Analytical analysis of the test steady-state thrust and power data results in 
functional relationships with blade angle. Taking derivatives with respect to 
blade angle provides the thrust and power derivatives. These derivatives 
could have been obtained graphically from curves of steady-state thrust and power 
against blade angle. These thrust derivatives are J = 0 values, but not the 
1 inear-with-J values. The tests were made at J = 0, thus the thrust starred 
value of equation (B18) remains unknown. For this case an estimate of the 
thrust starred derivatives is given in equations (B19) and (B20), in terms of 
thrust derivatives at J = 0. With the derivatives known and with equations 
(53) and (57), the cyclic-control values are determined from equations (73) 
through (76). Details of the theory application are given in the section titled. 
Comparison of Theory with Experiment. In that section, theory values of thrust 
change and pitching moment are compared with experiment, for various values of 
blade angle. The theory and test results correlate within the accuracy of the 
test measurements. Similar correlation is shown for the side force prediction. 
The theoretical ideal change of power is of the order 30 percent less than the 
total measured change of power. Extra power values are estimated in the section 
titled. Total Power Comparisons of Theory with Experiment, which shows good 
correlation between theory and experiment of total power change due to cyclic- 
control. In comparison with the four extra power estimates, the slipstream 
vibratory swirl power due to cyclic control appears as the largest. 


Isolated Cyclic-Control Forces and Moments 

In an aircraft control system it is desirable to initiate a given force or 
moment without other accompanying forces or moments. The cyclic-control counter- 
rotating propeller has this unique quality as shown in the development given in 
the section, Counterrotating Propeller. Qualitatively, isolation of forces and 
moments are demonstrated at the beginning of that section, and analytically in 
equations (124) through (127). The control is by means of the azimuth angle of 
each rotor. When ^02 of rotor 2 is set equal to ^01 of rotor 1, then the forces 
are zero and the moments are isolated. When T()2 is set equal to ?r + foi > then 
the moments are zero and the forces are isolated. 


Concurrent control of arbitrary force and moment of a counterrotating 
propeller . - Equations (118) through (121) apply for equal steady-state opera- 
tion of the clockwise rotation rotor 1, and the counterclockwise rotation 
rotor 2. The resultants of arbitrary moment and of force are given by 


C mn = < C m 2 + C n 2 >' /2 ' at t mn ’ 

where k " = 0 ' TV> 


k n 7r + tan 



(137) 
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2 2 1/2 -1 
Cyjsj ~ (Cy , dt ~ ky'ir + tan ^ 

c c 

uhere k Y . ^ (1 - -p^) 

The desired moment, force, and directions are C mn , Cy N , and fy N . The 

clockwise rotation rotor and the counterclockwise rotation rotor respectfully 
have the cyclic-blade angle variation given by 


+ Y-j COS (¥ - T Q1 ) 

(139) 

+ y 2 cos (¥ - Tq2 ) 

(140) 


For given values of Y-| , ¥q-j , Yp, and 4 , q 2> the C m , C n , Cy, and coefficients 
are determined from equations (118) through (121). With these values, the 
resultant moment and force are obtained from equations (137) and (138). In the 
application for isolating force and moment resultants from each other, the 
Y-angles become Y 2 = Yi, then equations (137) and (138) simplify to the solu- 
tions of equations (130) and (132). 

Multi propeller . - The forces and moments of a single propeller with cyclic- 
control cannot be isolated. In the section, Multi propellers In Line, it is 
shown how for two propellers in line as on a wing, the pitching moment and side 
and normal forces can be isolated, but yawing and rolling moments cannot. Three 
propellers in line is the minimum number for which all forces and moments can 
be isolated. For propellers in line the direction of propeller rotation of 
one propeller relative to the next, influences the magnitude of the isolated 
yawing and rolling moments. 


( 138 ) 


Effect of Cyclic Function 

Cyclic functions, f(T - t 0 ) , other than sinusoidal are investigated in the 
sections Cyclic function of cosine to m-power and Cyclic function for most 
moment with least extra power, also, with minimum blade pitching acceleration. 

The results show that when the cyclic function is a step function then the 
cyclic-control forces and moments are maximum, 27 percent larger than for the 
sinusoidal cyclic function. The ideal change in power is twice larger than that 
for the sinusoidal cycle. The step function means that blade angle is 3 0 + Y 
for half a cycle and 3o - Y for the second half, which results in an infinite 
blade angle acceleration at each change. The largest cyclic-control force or 
moment for a given ideal change of power is shown to be a cyclic function 
which is sinusoidal. Blade fatigue depends on the maximum blade angular accel- 
eration. The cyclic function which has the smallest maximum acceleration has 
ten percent less than that of the sinusoidal cycle and is given in equation (39). 
This function is very close to sinusoidal and the force or moment to power 
ratio is essentially the same as that of the sinusoidal cycle. 
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Propeller Inclination 


Propeller forces and moments and thrust and power changes are generated 
when the propeller is inclined to the free-stream velocity. These values are 
obtainable from the theory of reference 4. It is the derivatives with respect 
to blade angle of the steady-state thrust and power of the inclined propeller 
which are used for application in the present cyclic-control theory. The total 
propeller forces and moments are those due to cyclic-control plus those due to 
incl ination angle. 


CONCLUSIONS 


The equations for forces and moments, and thrust and power changes, due to 
a cycl ic-control led propeller, when analytically derived by power series, have 
led to a mathematical classic in simplicity. These cyclic-control equations 
are directly related to the derivative with respect to blade angle of the steady- 
state propeller thrust and power characteristics. The solution remains simple 
for any harmonic variation of cyclic-blade angle with azimuth angle during the 
blade revolution. An alternative sol ution, 1 imited to sinusoidal variation of 
cyclic-blade angle, is in terms of Bessel functions, which results in a closed- 
form solution for any value of maximum cyclic-blade angle. The study shows that 
a sinusoidal variation of cyclic-blade angle gives the maximum control force and 
moment for a given ideal change in power. A step variation of cyclic-blade 
angle gives the maximum control force and moment, which is 27 percent greater 
than that of sinusoidal variation, however, the ideal change in power is twice 
larger. A sinusoidal cycle with a small third harmonic is as effective as 
the sinusoidal alone, but has ten percent less maximum angular acceleration 
which may be a fatigue consideration. Analyses for estimating four extra power 
required terms show that the slipstream vibratory swirl excited by cyclic actions 
adds 15 to 20 percent to the ideal change in power, while several other contri- 
butors to power remain small. An example application of the theory to a cyclic- 
control test propeller shows good correlation between theory and test results of 
cyclic-control force and moment, and power and thrust change. 

Cyclic-control analysis is developed for the counterrotating propeller. It 
is shown that uniquely, a single counterrotating propeller can isolate all 
cyclic-control forces and moments, thus this propeller makes an ideal powered 
control, active control, and propulsive unit for a vehicle. A minimum of three 
and preferably four propellers in line, as on a wing, are needed to isolate all 
cyclic-control forces and moments. 

KENTRON INTERNATIONAL, INC. 

Hampton Technical Center 
an LTV company 
Hampton, Virginia 23666 
May 2, 1979 
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APPENDIX A 


QUASI -STEADY AZIMUTH ANGLE CHANGE DUE TO BLADE 
CYCLIC PITCHING AT CHORD POSITION f 

d 


For a blade that pitches about an axis x a distance from the leading edge, 
a linear with x velocity is induced along the blade chord, as shown in the 
drawing. This induced velocity creates 

blade 

IA—Il 

x a “ 

> V 


“a 


Ajixis of rotation 


"f — * 


T T ^ 


an induced camber on the blade which is 
positive for plus w a , and is a negative 
camber when a>a is negative. This camber N L 
changes the blade loading and lift, thus 


d£ 

dx 


_ 1 


p/b 


(<*r 




m 


The quasi-steady solution assumes that the lift is that with wa as a constant. 
The induced velocity is 



(A2) 


where b is blade chord, £ = x/b, £ a = x a /b. The thin-airfoil solution for 
lift of the blade with this induced velocity is 


C I = 2 " [ “o 




(A3) 


The value within the brackets represents the effective blade angle which for 
cycl ic control becomes 


3 0 + Yf(T - f 0 ) = 3 0 + Yf(Y - + 




where 


(A4) 


V 


R 


2imr 
COScJ) 9 


ph 


+ v 


con 


(A5) 


The angular velocity of a cyclic-control blade is the derivative of equation 
(A4) with respect to time, then 


_ dB o , dY , , df < 4 ' ’ V , b ,3 r . d “a 

“a ‘ dt" + dt f(T ‘ V + * at + ( T ’ 5 a> dt~ 
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which since g 0 , Y, and w a are assumed constant with time, becomes 


df(4< - 4\) d(4> - 4 , -| ) d-f ( v - 4 O 



Also, at r = R 

w ~ 3 r,r-| x ( \ 2 -| 1 / 2 l _ 3 n o 

V R = 2 wnR[l + (3^-) ] , b - j ttR g 

Combining equations (A6) and (A7) results in 


(A6) 


(A7) 


2tt | Y df(v - 4^) 


[i ♦ <£) 2 ] ,/2 


, 4* - 4S 


Equation (A8) is inserted into the effective blade angle equation (A4), then 


f(Y - O = ~ v-i) + 


t a ) df(T- T,) 
41 2 1/2 d( 4* - 4'-. ) 

[i + (f^) r 2 1 


which shows angular change due to the derivative term. For sinusoidal blade 
variation, f(f - 4 , -j) = cos(4' - 4'-|), then 

2 tt ~ (1 - E, ) 

COS ( 4* - 4' rt ) = COS(4' - 4',) sin(4' - 4',) 

+ (^ 2 1 1/2 


<! - ^ 
1 ♦ <s > 2 


C0S<¥ - + tan 


-1 2 » b (4 ~ e a ) 
[i + (^) 2 ] 1/2 


Since (o/Br is small, the factor of the cosine term is approximately unity, 
then with equations (1) and (A5), equation (A10) simplifies to 


v ax U a ) = -tan 


-1 2l B ^4 ~ 

D ♦ (^) 2 ] ,/2 
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The effective is Y ax (e a ) - f ax (|) > then 


ax 


= tan 


-1 


2 " f U a ’ 1 } 


[1 + 


/4J.\ 2-i 1 /2 

' 3-nr ' J 


{ AT 2) 


where £ a = x a /b is the blade chord position from blade leading edge about 
which the blade pitches, typically £ a = 1/2. The effective phase angle change, 
¥ ax , due to blade pitching at arbitrary values of £ a , will not be large since 
the ratio of propeller solidity to number of blades is small. 
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APPENDIX B 


APPROXIMATE STEADY-STATE PROPELLER THRUST AND 
POWER COEFFICIENTS AND DERIVATIVES 


From reference 4 the linear with advance ratio approximations of thrust 
and power coefficients are 


Cy* Gy( JqJ “ J)C0S( (3q - c) 

(Bl) 

C p “ e p( Jqp “ ^ ) sin 3^ 

( B2) 


where 


4.6 Ka 

e T Q 

1 1 + (3 + ~)ko 

(B3) 

10.4 Ka 
p i + a + 

(B4) 

J 0T = gtan 6 0 

(B5) 

g = 2. 1 + .5 tan 3| 

(B6) 

J 0P = J 0T + h 

(B7) 

8 d [1 + 7 sin 2 3(6 Q - 3-,)] 

^ 3 

15 sin 3 o cos 3 0 

(B8) 

d o " 6 o + - 8(6 o + B )sin B 1 

(B9) 


where k is blade section lift coefficient divided by 2ir, B is number of blades, 
6 0 is blade section drag coefficient at zero lift, and 3-j is given in equation 
(54). The nonlinear variation with advance ratio is approximated by a square 
root factor, as 


C 


T 



(BIO) 
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Cp* 

C P = 7172 

L p 

where 

e t - 1 + f e T 2 (' + I s 1" 2 s„) 2 0 - 4 )2 

E p - 1 + 9 e p 2 (l - /-) 2 
K p J 0P 


(Bll) 


(B12) 

(B13) 


Thrust Coefficient Derivatives 

The ratios of Cy'/Cy, and Cy /Cy are needed. The single prime and double 
prime indicate the first partial derivative and the second partial derivative 
with respect to g 0 . From equation (B5) 


OT 


OT 


2J 0T _ g 
sin2f3 o cos 2 e 


2J 


OT 


0 

2gsing 


0 


cos 2 g. 


3„ 

cos B. 


J o 0 

then the first derivative of equation (B1 ) , divided by Cy, is 


(B14) 

(B15) 


c f 


l (J ot - 0 J)s1n26 o - tan(S o - E > 


1/2 


The second derivative of equation ( B1 ) , divided by Cy, is 


(B16) 


C* r 2J nT sine 

_i = _i + yj — 

T L ( J q T " J)sin3 0 cos B Q cos(B Q - e). 

When advance ratio is zero, the advance ratio at zero thrust terms becomes unity 
in equations (B16) and (B17). 

Generally steady-state propeller thrust coefficients are presented in 
propeller charts as a function of advance ratio, J. Then the linear with 
advance ratio thrust coefficient is obtainable as 


1/2 


(B17) 
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(J 0T - J> 


(B18) 


8C t 

C T = 


J=J 


OT 


from which the derivatives with respect to s 0 or B can be evaluated. However, 
if only static conditions, J=0, data is available, then 3 Cj/ 3 J cannot be form- 
ulated. For this case the Cf derivatives can be obtained from equation (BIO). 
The first and second derivatives of C-f in terms of Cj, are 


2_= ( i[_ + iL) E 1/2 

c T c T e t )e t 


(B19) 


C f 


II II 

r C T E T C T 

[ r E T C T 


1(^L)2 + i!l_] E 1/2 

4 Ej ; 2 E t Jt T 


(B20) 


where Ey is given in equation (B12), and the derivatives of Ey are, for J = 0 


E T , = f e T 2(cos2e o + ll 81 " 4 ^) 


(B21 ) 


E T " = T e T 2( f cos40 o " sin 2 e o ) {B22) 

Application of equations (B19) and (B20) provides approximate predictions of 
the linear with J thrust coefficient derivatives in terms of steady-state 
test or propeller theory thrust coefficient derivatives. 


Power Coefficient Derivatives 

The derivative ratios Cp'/Cp, and Cp''/Cp are obtained by taking the deri- 
vatives of equation (Bll) with respect to Bo* The first and second derivatives 
divided by Cp are 


C P 

T 

ifp 

(B23) 

C P 

C P 

2 E p 

Cp' 

Cn 

C** 1 

cs 

E P C p' 3, E F\ 2 1 E p' 
" En c* + ^lj'2 E n 

(B24) 
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where E p is given in equation (B13), and the C* and E p derivatives are 


c r 

c ? 


i + j op 


tan$ 


J 0P ' J 


( B25 ) 


C P* ' = •, + J 0P 2J 0P 

C P " J 0P ■ J ( J 0P ■ tan ^o 


(B26) 


EA = 18 e n 2 (1 - ^~) ^ 

J 0P J 0P 


■P J 


E p ' = 18 e 


OP 
2 J 


J \ J 0P 


p D-C(l (2-3j^-)(^) 2 ] 

r u nn u An n U o p ^QP 


'OP U 0P U 0P 

where Jq P is given in equation (B7), and the derivatives are 


JAn = 


2J 


OT 8d o J 21 sin6 ^ 0 ' 6 1^ 1-4 sin2e ( 

+ 


'OP sin2g 15 j . 3 

p o [ sine Q cos 3 0 

[1+7 sin 2 3(3 q - 3^]! 


. 2 4 

sin 3 0 cos 3 0 


,, 2J 0T 16d o f 63 cos6 (e o " 3]) 21 (1 - 4 sin 2 3 Q ) 

j op " 2„ + 15 ) ~7 777 77 X 7“4 sin6 (3 0 - b-j) + 


(B27) 


(B28) 


(B29) 


cos 3 


o 


sin3 0 cos 3 q 


sin 3 0 cos 3 0 


1 - 3 sin 2 3 Q + 8 sinV 


- [1 + 7 sin 2 3(3 q - e n )]' 


. 3 5 

sin 3 0 cos 3 0 


(B30) 


At zero velocity, or in hover, the advance ratio is zero, then E p = Ep 1 = 0, 
and the power coefficient derivatives greatly simplify. For J = K 0 


r 1 ii 

_P = _J + ^OP 

C P " tane o J 0P 



+ 2J 0P + ^OP 
J 0 p tans 0 J op 


(B31) 
( B32) 
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